Critical phenomena in networks 
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We develop a phenomenological theory of critical phenomena in networks with an arbitrary 
distribution of connections P{k). The theory shows that the critical behavior depends in a crucial 
way on the form of P{k) and differs strongly from the standard mean-field behavior. The critical 
behavior observed in various networks is analyzed and found to be in agreement with the theory. 
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Many complex, interacting systems such as the brain, 
Internet, social systems, etc. are recognized as networks. 
They demonstrate a spectrum of unique effects [1-5], 
which distinguish them from all other known structures 
of condensed matter. One could expect that such sys- 
tems have peculiar cooperative phenomena. Only the 
first studies of particular cooperative models on particu- 
lar networks were performed recently [6-18]. They were 
focused on traditional models in statistical physics, such 
as the Ising and XY models [7-14], percolation [15-17], 
epidemic spreading, et al. in scale-free networks [18]. 
It was revealed that their critical behavior is more rich 
and extremely far from that expected from the standard 
mean-field theory. A deviation from the usual mean-field 
critical behavior appears when, depending on the model, 
the fourth, (fc**), or third, {k^), moments of the degree 
distribution P{k) diverge. Here, degree is the number of 
connections of a vertex. When (fc^) diverges, all mod- 
els undergo an unusual phase transition of infinite order. 
This case is of primary importance because many real 
networks such as the Internet and biological nets are de- 
scribed by P(fc) with infinite (fc^). 

Why do critical phenomena in networks differ so much 
from those in usual substrates and what is their common 
origin? Why do all investigated models demonstrate uni- 
versal behavior when (fc^) diverges? In order to answer 
the questions raised above and analyze results of simula- 
tions and experiments from a general point of view, it is 
necessary to have a general theory which is not restricted 
by specific properties of any model. 

The advantage of the phenomenological approach is 
that the origin of interactions and nature of interacting 
objects are irrelevant. These can be spins, percolating 
clusters, biological objects and etc. In this Letter, we de- 
velop a phenomenological theory of cooperative phenom- 
ena in networks. Our approach is based on concepts of 
the Landau theory of continuous phase transitions. It is 
shown that the critical behavior in networks has a univer- 
sal character and is determined by two general properties: 
(i) the structure of a network and (ii) the symmetry un- 
derlying a model. We find that in networks described by a 



degree distribution with divergent moments, the thermo- 
dynamic potential of an interacting system with cooper- 
ative effects is a singular function of an order parameter. 
The theory is in complete agreement with exact results 
for the Ising and Potts models and percolation on 'equi- 
librium' uncorrelated random networks. On the basis of 
the theory we discuss results of theoretical studies and 
simulations of the Ising and XY models and epidemic 
spreading on evolving and small- world networks. 

Let us consider a system of interacting objects. Inter- 
actions or links between these objects form a net. We 
assume that some kind of 'order' can emerge. This 'or- 
dered' phase may be characterized by some quantitative 
characteristic x while it will vanish in a 'disordered' phase 
above a critical point. In order to study the critical be- 
havior we assume that the thermodynamic potential $ of 
the system is not only a function of the order parameter 
X but also depends on the degree distribution: 

oo 

<^>{x,H) = -Hx + ^P{k)(f>{x,kx) . (1) 
fc 

Here H is a field conjugated with x. It should be empha- 
sized that Eq. (1) is not obvious a priori. The function 
4>{x, kx) may be considered as a contribution of vertices 
with k connections. There are arguments in favor of this 
assumption. Let us consider the interaction of an arbi- 
trary vertex with k neighboring vertices. In the frame- 
work of a mean-field approach, k neighboring sites having 
a spontaneous 'moment' x produce an effective field kx 
acting on the site 0. This indicates that the expansion is 
over not only x but also kx. 

We assume that 4'{x, y) is a smooth function of x and 
y and can be represented as a series in powers of both x 
and y: 

oc oo 

^^^^^xV, (2) 

m=0 (=0 

where ipmi are functions of 'temperature' T and H . We 
use the term 'temperature' for convenience. If a change 
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of another parameter leads to the emergence of 'order', 
then this control parameter can be used as T. 

There are general restrictions on the coefficients cpmi- 
For an arbitrary P(fc), at zero field H = 0, the expansion 
of <^{x,H = 0) over x must contain only the second and 
higher powers of x. Therefore, 0oo = 4>Qi = 4'ia = 0. 

$ must be finite if (k) is finite. This condition is sat- 
isfied if at y 1 and an arbitrary x the function (p^x, y) 
increases slower than y: 

(t){x,y) ^ g{x)y , (3) 

where g{x) is a function of x. 

Eqs. (1) (3) are the basis of our theory. 

In the framework of the Landau theory, it is assumed 
that the thermodynamic potential $ near a phase tran- 
sition point can be represented as a series in powers of 
the order parameter x: 

^{x, H) = -Hx + f2x'' + hx^ + Ux^ + .... (4) 

Higher terms in the expansion (4) are irrelevant for the 
critical behavior. The coefficients /„ are related to the 
coefficients (pmi in the expansion (2). Let us suppose first 
that the critical temperature is finite. (The case of 
an 'infinite' critical temperature will be considered sep- 
arately.) Following Landau, we assume that near the 
critical temperature the coefficient /2 can be written as 
a(T — Tc) where a > due to the condition of stability of 
the disordered phase at T > Tc- /2 = at the transition 
point. This gives the following equation for Tc- 

<p2o{Tc) + <Pn{Tc) {k) + cj>02{Tc) (fc') = . (5) 

According to this equation, in the case ^20 = the criti- 
cal temperature Tc depends only on the ratio (fc^) / {k). 
Such a dependence was revealed in all studied models on 
networks [8,9,15,17,18]. In the ordered phase the order 
parameter x is determined by the condition that ^{x, H) 
is minimum; 

d^{x,H)/dx = Q. (6) 

Solving Eq. (6), we find x{T,H) and all other physi- 
cal parameters: the response function Xx = dx/dH, the 
specific heat C = -Td?^/dT'^, et. al. 

The condition of the stability of the ordered phase 
near the continuous phase transition demands that ei- 
ther (/a) > or if (/g) = 0, then {U) > 0. In 
the general case, /2n+i 7^ 0. If for symmetry reasons 
^{x,H) = -iJ), then /2n+i = 0. This condition 

takes place if 0™/ = 0atm-|-l = 2n-|-l. 

Let us consider analytical properties of $ in the gen- 
eral case. The n-th derivative of ^{x,H) at a; = is 
equal to 

n 

$(")(0) = d^<^{x,H)/dx"l^, = n\J2<t>n-i,i {k') . (7) 



If all the moments ^fc') = J2k -^i^)^^ converge, then all 
the derivatives $("^(0) are finite. Therefore, $(x,ff) is 
a smooth function of a;, at least at small x. However, in 
most interesting real networks, P{k) has divergent mo- 
ments. This leads to a very important consequence. If 
moments (fc') with I < p converge and moments (fc') 
with I ^ p diverge, then $^'^(0) with I < p arc finite 
while $(')(0) 00 ioT I ^ p. In this case ^{x,H) is a 
singular function of x and can be written as follows: 

p-i 

^{x,H) = -Hx + ^fn{k)x''+xPs{x), (8) 

where s{x) is a singular function: s(0), s*^")(0) — > cxd, but 
a;s (a;) —> at a; —> 0. This result is of primary impor- 
tance. It is the singular function that can lead to a devi- 
ation from the standard mean-field behavior. A method 
for determining s{x) is proposed below. For convenience 
we use a power-law degree distribution P(k) oc k^^. 
Then, (fc'*) diverges for 7 < 5, (k^) diverges for 7 < 4, 
and (A;^) is divergent for 7 < 3. 

We begin with the case /2n-i-i = 0. Let us consider the 
critical behavior for /4 > and diff'erent 7. 

(a) 7 > 5. Here, (fc^) converges. The coefficient /4 is 
finite and $ has a usual form (4). At H = Eq. (6) 
leads to the standard critical behavior: 

a; ocT^/^ AC ^ 0, xa;(i? = 0) oc (9) 

Here, t = 1 - T/T^, AC - 1/ (/c^) is the jump of the 
specific heat at Tc- It tends to zero when 7-^5. 

(b) 3 < 7 ^ 5. Here, {k'^) diverges. According to Eq. 
(8), a singular term a;*s(a;) is expected to appear in 
The function (f){x, y) in Eq. (2) may be divided into two 
parts: 

00 

(l>{x,y)=Y, <l>mix^y' + G{x,y). (10) 

m=0 l^i<7-l 

The first term determines a non-singular contribution to 
$ in Eq. (1). The second term, G{x,y), gives the sin- 
gular contribution. For 4 < 7 ^ 5 the function G{x,y) 
has the following properties: G(0, y) oc at small y, 
(i"G(0,y)/dy" 9^ for n > 4 at 2/'= 0. Owing to the 
condition (3), the function G{x,y) increases slower than 
j/^ at 2/ > 1. 

It is convenient to use a continuous approximation for 
a degree distribution P{k) = Ak~'^, where A is a nor- 
malization factor. Then Eq. (1) may be rewritten in the 
following form: 

f°° dk 

^{x,0)=hx^ + ipix)x''+A G(x,xk), (11) 

Jm 

where m is the smallest degree in P{k), ip{x) is a smooth 
function which is determined by the convergent moments 
(fc') with I ^ 3, ip{Q) ^ 0. In order to find $ at small x, 
one can put G(a;, y) w G(0, y) into Eq. (11). Considering 
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the integral over a variable y = xk, one can show that 
the region mx < y gives a leading contribution, b is 
a model parameter. We obtain 

^x, H) = -Hx + a{T - Tc)a;2 + Ax'^ ln{b/x) , (12) 

^x, H) = -Hx + a{T - T^)x'^ + Bx'^'^ , (13) 

for 7 = 5 and 3 < 7 < 5, respectively. A and B are 
constants at T = Tc. The critical behavior for the case 
/a = 0, /4 > at 7 > 3 is summarized in Table I. It was 
observed in the Ising model on 'equilibrium' uncorrelated 
random networks [8,9]. 

Note that the divergence of {k'^) does not change the 
critical behavior of Xx- 

(c) 2 < 7 ^ 3. Now (A;^) diverges. Equation (5) shows 
that in this situation there is no phase transition at any 
finite temperature. According to Eq. (8), a singular term 
x'^s{x) is expected to appear in $. Calculations of $ may 
be carried out in a similar way as above, using the func- 
tion G{x,y) = ^{x,y) - Em^i^rnix^'v, Eq. (10). We 
find 

<P{x,H) = -Hx + Cx"^ - Dx^lii{b/x) , (14) 

^{x, H) = -Hx + CV - D'.x'^-i (15) 

for 7 = 3 and 2 < 7 < 3, respectively. The coefficients 
C,D,C' and D' are functions of T and H. If C, D , C 
and D' arc positive, then at arbitrary finite temperature 
T 1 in zero field H = 0, Eq. (6) has stable nontrivial 
solutions: 

X = 6exp(-(2(7 + D)/{2D)), Xx = l/(2£') , (16) 

V(7-l)£'7 ' = 2(3 - 7)C' ^^^^ 

for 7 = 3 and 2 < 7 < 3, respectively. As in the Landau 
theory, the results (16) and (17) are obtained in terms 
of the coefficients of the thermodynamic potential. Note 
that in the situation where a phase transition is absent 
at any finite temperature, the phenomenological theory 
can not determine the temperature behavior of the coef- 
ficients C, -D, C" and D' . In this situation, the temper- 
ature dependences of the coefficients can be found only 
by a microscopic theory (see below). 

In the general case the symmetry of the model admits 
non-zero . In the case /a > the analysis of the analyti- 
cal properties of <f> can be performed as above. For 7 > 4, 
$ is a smooth function of x and leads to the standard 
critical behavior. At 7 ^ 4, $ contains a singular term: 
(a) In a; at 7 = 4; (b) x'^~^ in the ranges 4 > 7 > 3 
and 3 > 7 > 2; (c) a;^ Ina; at 7 = 3. The corresponding 
critical behavior for 7 > 3 is represented in Table I. This 
behavior was observed for percolation on uncorrelated 
random networks [17]. In the range 2 < 7 ^ 3 when (A:^) 
diverges, $ has the universal form, Eqs. (14) and (15). 

In the case /a < 0, or < if fz = 0, when (A:^) 
converges, a first-order phase transition occurs. In agree- 
ment with this prediction we found such a transition in 



the (7— state Potts model with (7 ^ 3 on 'equilibrium' un- 
correlated random networks by use of the approach of 
Ref. [8]. In the limit 7 — > 3 the jump of the order param- 
eter at the transition tends to zero, and the transition 
transforms into the infinite order phase transition. A 
detailed study of the transition will be given elsewhere. 

In order to complete Table I, let us discuss the temper- 
ature behavior in the case 2 < 7 ^ 3 within the micro- 
scopic theory of the Ising model and percolation on 'equi- 
librium' uncorrelated random networks [7-9,17]. For this 
purpose we use the more general ferromagnetic q'— state 
Potts model which aX q = 1 and 2 is equivalent to perco- 
lation and the Ising model, respectively (see, for example 
Ref. [19]). Using the approach of Ref. [8] we obtain that 
in the Potts model a continuous phase transition occurs 
at the exact critical temperature 

^^-2/ In ^p^_2(fc) • (18) 

Hereafter, we set the energy of ferromagnetic interac- 
tion between nearest neighbors J — 1. The parameter 
Pc — 1 — exp(— 2/Tc) determines the percolation thresh- 
old [16,17] and establishes the relation between the Potts 
model and percolation. 

Let us consider the case 2 < 7 ^ 3. Here, (A^) diverges. 
Tc is infinite for the infinite networks. In any finite net- 
work, (fc^) < 00, and Tc is finite, although it may be 
very high, Tc = 2{k'^) / {{k) q). In the temperature re- 
gion T ^ 1, where a; ^ 1, the Potts model demonstrates 
universal behavior at all > 1: 

x«(g/(fc))e-«^/W, XxOcT-i, (19) 
X ~ r-i/(3-^), Xx oc T-2 (20) 

for 7 = 3 and 2 < 7 < 3, respectively. Without the 

continuum approximation, we have instead of (fc) in the 
exponential, a constant which is determined by the com- 
plete form of P{k). In Ehrenfest's classification, this 
transition is of infinite order, as all temperature deriva- 
tives of $ are finite. The results (16) and (17) agree with 
Eqs. (19) and (20) if we put C, C" oc T^, D, D' oc T. The 
exponential behavior (19) has been revealed in epidemic 
spreading within scale-free networks with 7 = 3 [18] and 
in percolation on these networks [17]. 

At small X there is the following relationship be- 
tween the response function Xx and the susceptibility 
X = dM/dH: x ~ 2/gT + (fc) Xx/q- At 7 > 3 the crit- 
ical behavior of x is determined by Xx- X ~ {k)Xx/<l- 
At 2 < 7 ^ 3 we have x oc l/qT, as the paramagnetic 
contribution 2/qT is of the order of Xx or larger. 

Let us discuss the results of theoretical studies and 
simulations of critical phenomena in different networks 
on the basis of the phenomenological theory. 

The phenomenological theory as well as the Landau 
theory assumes that the contribution of fluctuations to 
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the thermodynamic potential is small. Above we have 
shown that the theory gives the exact critical behavior 
of the Ising model, percolation and Potts model on 'equi- 
librium' uncorrelated random networks [8,9,17]. The rea- 
son for this is that these networks have a local tree-like 
structure and vertices are statistically equivalent [20-22] . 
Due to these properties, vertices in the networks can be 
regarded as forming a Bethe-lattice structure for which a 
mean- field approach is exact [23] . It means that the fluc- 
tuation contribution is negligibly small. Note that in a 
graph with a Cayley tree-like structure, vertices are sta- 
tistically inequivalent, and a mean- field approach is valid 
only for vertices deep within such a graph [23] . 

The origin of a deviation from the standard mean-field 
critical behavior is different for regular lattices and net- 
works. In a regular lattice the deviation is caused by 
strong critical fluctuations which depend crucially on the 
space dimension. In networks the deviation is brought 
about by the most connected vertices which induce strong 
correlations in their close neighborhood. With decreas- 
ing exponent 7 the relative number of highly connected 
vertices increases and their role turns out to be more 
important. 

Recent simulations [7] of the Ising model on the 
Barabasi- Albert scale- free network with the degree distri- 
bution P{k) = Ak~^ revealed a temperature behavior de- 
scribed by Eq. (19) as well as in the Potts model on 'equi- 
librium' uncorrelated random networks. Unlike the latter 
networks, the Barabasi-Albert net is correlated. Never- 
theless, results of the simulations agree with the universal 
temperature behavior predicted by the phenomenological 
theory for networks with the divergent moment (fc^). 

The exact results for percolation on small-world net- 
works were obtained in Ref . [6] . The Ising and XY mod- 
els on small- world networks were studied analytically and 
by use of Monte Carlo simulations [10-13]. These net- 
works, introduced by Watts and Strogatz [24], have a 
Poisson-like degree distribution and large clustering co- 
efficients. It was found that the critical behavior is char- 
acterized by the standard mean-field critical exponents. 
This result is in agreement with the prediction of the phe- 
nomenological theory that the standard mean-field crit- 
ical behavior should occur in networks described by a 
degree distribution with convergent moments. 

Thus, the analysis of critical behavior of the Ising, 
Potts and XY models, percolation and epidemic spread- 
ing on uncorrelated random, scale-free and small-word 
networks shows that the critical behavior in networks de- 
pends crucially on the form of the degree distribution and 
the symmetry underlying a model in a complete agree- 
ment with the phenomenological theory. The studied 
networks differ in clustering coefficients, degree correla- 
tions, etc. However, the situation looks like these char- 
acteristics are not relevant for critical behavior. Their 
role in the formation of the critical exponents is still un- 
clear. Further investigations in this topic are necessary. 
It would be interesting to find a network where criti- 
cal behavior differs from a mean-field one due to strong 



fluctuations. However, even in this case it is expected 
that when T tends to Tg, a temperature region of the 
mean-field critical behavior precedes a region of strong 
fluctuations. 

Real systems like the Internet, WWW or biological 
networks have a network structure with a power-law de- 
gree distribution with exponent 7 below 3, see Refs. [3,4]. 
The phenomenological theory predicts for this case the 
power-law critical behavior (20), see also the last line in 
Table I. This behavior agrees with an empirical obser- 
vation [5] for the nd.edu domain of the WWW, where 
the variations of the size of the giant component under 
random damage were studied. This size and the number 
of damaged vertices play a role of the order parameter x 
and the control parameter, respectively. 

In conclusion, the phenomenological theory shows that 
the deviation of the critical behavior of interacting sys- 
tems with a network structure from the standard mean- 
fleld behavior emerges when a degree distribution has di- 
vergent moments. The theory predicts different classes of 
critical behavior in agreement with microscopic studies of 
the Ising and Potts models, and percolation on 'equilib- 
rium' uncorrelated random networks. It also agrees with 
results previously obtained for various models on small- 
world and evolving networks. The theory can easily be 
generalized for models with a multicomponent order pa- 
rameter and can give useful insight into collective effects 
in different networks discussed in connection with the 
Internet, biological networks, etc. Using this approach, 
one can also study the critical relaxation in models on 
networks. 
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5C{T < Tc) 



/3=0, 
/4>0 



/3 >0 



arbitrary 
fs and /4 



7 > 5 

7 = 5 
3 < 7 < 5 
7 > 4 
7 = 4 
3 < 7 < 4 

7 = 3 
2 < 7 < 3 



^1/2 

ri/V(lnr-i)i/2 

.^1/(7-3) 



r/(lnr-i) 

^1/(7-3) 



j.-l/(3-7) 



jump at Tc 
1/lnr-i 

7-(5-7)/(7-3) 
T 

T/(lnV~i) 

^(6-7)/(7-3) 
J-(7-l)/(3-7) 



TABLE I. Critical behavior of the order parameter x, the specific heat SC, and the response function of a model on 
networks with a degree distribution P{k) ~ k ' for various values of the coefficients /s and fi, and exponent 7. t = 1 — T/Tc, 
c is a constant which is determined by the complete form of P{k). In the case fs < 0, or /4 < if /s = 0, at 7 > 3, the system 
undergoes a first-order phase transition. 
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